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ABSTRACT 



The previously published atmospheric neutrino data did not distinguish whether muon 
neutrinos were oscillating into tau neutrinos or sterile neutrinos, as both hypotheses fit 
the data. Using data recorded in 1100 live-days of the Super-Kamiokande detector, we 
use three complementary data samples to study the difference in zenith angle distribution 
due to neutral currents and matter effects. We find no evidence favoring sterile neutrinos, 
and reject the hypothesis at the 99% confidence level. On the other hand, we find that 
oscillation between muon and tau neutrinos suffices to explain all the results in hand. 



I. INTRODUCTION 

In previous papers pHq] Super-Kamiokande reported evidence for the oscillation of muon neutrinos 
produced in cosmic-ray induced showers in the atmosphere. This evidence rests largely upon a strong 
zenith angle dependent deficit in the muon data, which does not appear in the electron data, and hence 
limits the amount of oscillation into electron neutrinos. In fact, the results only demonstrate that muon 
neutrinos are disappearing, depending on the energy and flight distance. The most plausible scenario is 
that muon neutrinos oscillate to tau neutrinos, most of which are below the 3.4 GeV neutrino energy 
threshold for charged current tau production. The few charged current tau events created (we expect 
approximately 65 in our current sample) typically fail our cuts that identify a single electron or muon; 
indeed, due to the high energy of the interaction and the multiple decay modes of the tau, it is difficult 
to isolate any set of events that can be uniquely identified as due to charged current v T interactions. 

An alternative scenario that can explain muon neutrino disappearance is oscillation with a sterile 
neutrino iy s ), so named because it has neither charged current (CC) nor neutral current (NC) inter- 
actions. Neutrino oscillation has also been employed to explain two other experimental anomalies: the 
long-standing deficit of solar neutrinos g, and the appearance of electron antineutrinos in the LSND 
experiment ||. The three oscillation signatures, LSND, atmospheric, and solar, are manifested by three 
widely separated values of the mass-squared difference, Am 2 = mf — m 2 . Because Am^ must equal 
Am 2 2 + Am 2 3 , au three signatures cannot be accommodated with three neutrino states. Any additional 
light neutrino must be sterile to satisfy the well-known bound of three neutrino flavors that couple to the 
Z°§. _ 

In this letter, we use more than 1000 days exposure of atmospheric neutrino data collected by the 
Super-Kamiokande detector to distinguish the behavior of v^ <-» v r oscillation from v^ <-» v s oscillation. 
First, for v^ <-> v s oscillation, one should observe fewer neutral current events than for v^ <-» v r oscillation; 
in the latter case the neutral current event rate is unchanged by oscillations. 

Second, the interaction of the neutrinos with matter leads to a difference in the oscillation probability. 
The coherent forward scattering of v r and v^ are identical, therefore the presence of matter in the neutrino 
path does not modify the oscillation probability for a neutrino of energy E v that travels a distance L in 
vacuum: 

P v ^ Vr = sin 2 20„sm 2 (7r~), (1) 

where sin 29 v is the mixing angle between the two neutrino states, and I in vacuum is given by l v = 
4irE„/ Am 2 . In contrast, v s does not interact with matter by definition, whereas v^ does interact with 
matter via the neutral current. This introduces an effective potential which modifies the mixing angle 
and oscillation length J7J : 



sin 2 °™ = 77 — .._„„ ^ , _-__2 nn > ( 2 ) 



sin 2 29 v 
(C - cos20„) 2 + sin 2 29 v ' 

/,„ = ^ (3) 

(C- cos 20„) 2 + sin 2 29 v 



The parameter £ is given by ^fy/^E^GFNn/Am 2 , where N n is the neutron density in the matter traversed 
by the neutrino, the minus sign is for neutrinos, and the plus sign is for antineutrinos. For the density 
of matter in the earth, £ reaches unity for E v of 5 GeV x Am 2 /(10 _3 eV 2 ). The Super-Kamiokande data 
indicates a likely value for Am 2 of 3 x 10~ 3 eV 2 , which means that neutrinos with energy greater than 
approximately 15 GeV will have the oscillation probability suppressed by matter effects if the oscillation 
is Vp <-> v a . Neutrinos of lower energy will have approximately the same oscillation probability in matter 
as in vacuum, even if the oscillation is v^ «-> v s . 

II. THE DATA SETS 

Super-Kamiokande is a 50 kilotons water Cherenkov detector employing 11,146 photomultiplier tubes 
(PMTs) to monitor an internal detector (ID) fiducial volume of 22.5 kilotons. Incoming and outgoing 
charged particles are identified by 1885 PMTs in an optically isolated outer volume (OD). Details of the 
detector, calibrations, and data reduction can be found in Refs. |y-§|. Super-Kamiokande has collected 
9178 fully-contained (FC) events and 665 partially-contained (PC) events in a 70.5 kiloton-year (1144 
days) exposure, starting in April 1996. FC events deposit all of their Cherenkov light in the ID while 
PC events have exiting tracks which deposit some Cherenkov light in the OD. The vertex position 
is reconstructed and the number of Cherenkov rings are counted using PMT pulse height and timing 
information. The directions and momenta are reconstructed and the particle types are identified as "e- 
like" or "//-like" for each Cherenkov ring. In the current FC sample, there are 3107 single-ring e-like 
events, 2988 single-ring /i-like events and 3083 multi-ring events. 

This detector has also collected 1269 upward through-going muon (UTM) events during 1138 live 
days. Muons which leave both entrance and exit signal clusters in the OD are regarded as through-going 
muon events; those that are upward going are produced by atmospheric neutrino interactions in the 
surrounding rock. We required a minimum track length of 7 meters in the inner detector (corresponding 
to a minimum muon energy of 1.6 GeV) and a zenith angle cos© < (cosO = — 1 means vertically 
upward-going events). Because of finite fitter resolution and multiple Coulomb scattering of muons in 
the nearby rock, some down-going cosmic-ray muons appear to be coming from —0.1 < cos© < 0. This 
background was estimated to be 9.1 ± 0.8 events M, which was subtracted from the most horizontal bin. 

III. ANALYSIS 

A. Fully Contained Single-Ring Data 

First, utilizing only the FC single-ring events, we have examined the hypotheses of two-flavor v^ «-► v T 
and v^ <-» v s oscillation models using a \ 2 comparison of our data and Monte Carlo (MC), allowing all 
important MC parameters to vary, weighted by their expected uncertainties. For v^ <-» i/ s , the effects of 
matter on neutrino propagation through the earth were taken into account by a numerical evolution where 
the density of the earth was divided into 94 discrete steps in radius based on Ref. |§ . Furthermore, matter 
effects are different with positive or negative Am 2 except for full mixing case. Therefore, we evaluate 
three models of oscillation, (a) v^ <-> v T , (b) v^ «-► i/ s (Am 2 > 0), and (c) v^ <-> i/ s (Am 2 < 0). The data 
were binned by particle type, momentum, and cosO. A \ 2 is defined as: 

2 _ ^ {N data -N MC (sin 2 26,Am 2 ,6 3 )\ 2 ^ ( e^ 2 

Xfc= ^ V ~o ) M^ 

where the sum is over five bins equally spaced in cos O and seven (six) momentum bins for e-like (/i-like) 
events. N^ata is the measured number of events in each bin, a is the statistical error, and Nmc is the 
weighted sum of MC events. The definition of x 2 : an d the treatment of systematic uncertainties, ej, 



is identical to that in Ref. Q, except we exclude the PC events (which we employ later in the present 
report). 

The best-fit values of oscillation parameters are summarized in Table Q. With the best fit parameters 
for v^ «-> v T , we expect only 16 single-ring events from v T charged current interactions in the current 
sample. Moreover, matter induced modifications to oscillations do not produce significant effects due to 
the relatively small energy (~ 1 GeV) of the parent neutrinos for the FC events. Therefore these three 
hypotheses for oscillations are indistinguishable by this data sample alone. 

TABLE I. Best fit oscillation parameters for fully contained sample. 



Mode 


Am 2 (eV 2 ) 


sin 2 26» 


xLin/d-O.f 


^ M <-» Vt 

u M <-» ^ s (Am 2 > 0) 
Vy. <-> u e (Am 2 < 0) 


3.2 x 10" 3 

4.0 x 10" 3 
3.2 x 10" 3 


1.000 
0.995 
1.000 


61.33/62 
62.56/62 
62.62/62 



B. Multi-Ring Sample 

Next we employ a neutral current (NC) enriched sample of events selected from multi-ring (MR) data. 
By definition a sterile neutrino does not interact with matter even through the neutral current channel, 
while a tau neutrino continues to experience the same neutral current interactions as did the original 
muon neutrino. We measure the zenith angle distribution of NC events to distinguish between v^ *-> v T 
and i/p <-> v s : if pure v^ <-> v r oscillations are operating, then the up/down ratio should be nearly unity; 
if i/ M <-> v s oscillations dominate, the up/down ratio will be measurably smaller. 

In order to obtain a sample with an enhanced with NC events, we applied the following selection 
criteria: (1) vertex within the fiducial volume and no exiting track; (2) multiple Chcrenkov rings; (3) 
particle identification of the brightest ring is e-like; (4) visible energy greater than 400 MeV. 

The first criterion provides a contained event sample, the second and third criteria serve to enrich the 
NC event fraction. The fourth criterion helps to obtain good angular correlation between the incident 
neutrino and the reconstructed direction, defined as the charge weighted sum of the ring directions. The 
mean angle difference between the parent neutrino and reconstructed directions is estimated to be 33°. 
According to our MC study, for no oscillations (and v^ *-> v T oscillations at best fit parameters), the 
resultant fraction of NC events is 29% (30%), u e CC is 46% (48%) and ^CC is 25% (19%) (and ^ T CC 
is 3%). In contrast the FC single ring sample contains only ~ 6% NC events. In the current exposure, 
1531 events satisfy the above criteria. Figure [y(a) shows the zenith angle distribution of these events 
with predictions from the MC. 

We utilize an up-to-down ratio as the discriminant, which cancels some systematic uncertainties (oth- 
erwise dominated by the large uncertainty in absolute rates). In this context we define "upward" as a 
cosine of zenith angle less than —0.4, and "downward" as greater than +0.4. There are 387 upward events 
and 404 downward events. Figure 0(b) shows the Am 2 dependence of the expected up/down ratio in the 
case of full mixing (sin 2 29 = 1). For Am 2 of 3.2 x 10 -3 eV 2 , the data are consistent with v^ «-> u T , while 
the data differ from the prediction for v^ <-> v s oscillation by 2.4 standard deviations. 

We estimated the total uncertainty in the up/down ratio of the data and MC to be ±2.9%, dominated 
by the 2.6% uncertainty in the neutrino flux caused by the absorption of muons in the mountain above the 
detector. All other sources of systematic uncertainty such as background contamination, the up/down 
response of the detector, the uncertainty of the NC cross sections, and the difference between the up/down 
ratio of two independent flux calculations pffOj are less than 1%. 

C. Partially Contained Sample 

Next we report on the search for matter effects by using high energy partially contained events. As 
discussed above, matter effects impact only v^ <-> v s oscillations, where at high energies the matter effect 
suppresses oscillations. Partially contained events in Super-Kamiokande are estimated to be 97% pure 
v^ charged current, with a mean neutrino energy of 10 GeV. In order to select higher energy v^ events, 
which are more sensitive to matter effects, we additionally require visible energy greater than 5 GeV. 
We estimate the typical energy of the parent atmospheric neutrino is 20 GeV. After cuts are made upon 
the current data sample we find 267 events. Figure [j](c) shows the zenith angle distribution of these 
events with predictions from MC, as before. Again we employ an up/down ratio to cancels systematic 
uncertainties, with the same angular definition as used for the multi-ring sample. There are 43 "upward" 
events and 84 "downward" events. Figure |l|(d) shows the Am 2 dependence of expected up/down ratio 
in the case of full mixing. For Am 2 of 3.2 x 10~ 3 eV 2 , the data are consistent with v^ <-> v T oscillation, 
whereas it differs from v^ «-v v s oscillation by 2.3 standard deviations. 

We estimated the total systematic uncertainty in the up/down ratio to be ±4.1%, dominated by the 
3.4% uncertainty caused by the mountain above the detector and the 2.0% uncertainty caused by possible 
background contamination by cosmic-ray muons. All other sources of uncertainty were less than 1%. 



D. Upward Through-going Muon Sample 

Next we report on the search for possible matter effects by using upward through-going muon events. 
The approach of this analysis is similar to that for the PC events. Because the typical energy of the UTM 
parent neutrino is approximately 100 GeV, matter effect suppression should appear most prominently 
in this data setjj. In the current sample, we have 1259.9 events after background subtraction. Figure 
0(e) shows the zenith angle distribution of these events with predictions. Again we utilize a ratio as the 
test parameter, dividing "vertical" and "horizontal" at cosine of zenith angle = —0.4. Figure ll](f ) shows 
the Am 2 dependence of the expected vertical/horizontal ratio in the case of full mixing. At the point of 
3.2 x 10~ 3 eV 2 , the data are consistent with ja, <-> v T oscillation, while i/„ <-> v s oscillation differs from 
the data by 2.9 standard deviations. 

We estimated the total systematic uncertainty in the horizontal/vertical flux ratio to be ±3.3%, domi- 
nated by the 3% uncertainty in the ir/K production ratio in the cosmic-ray interaction in the atmosphere. 
All other sources ofsystematic uncertainty, including the background contamination in the most horizon- 
tal bin |3|], th e sp ectral index of the neutrino flux, and the difference between two independent flux 
calculations |J9|,|l0f were 1% or less. 

E. Combined Analysis 

Finally, we performed a combined statistical analysis of the multi-ring, high energy partially contained, 
and upward-going muon data sets. For each sample (i = MR, PC, UTM) we construct a one degree of 
freedom \ 2 defined by: 

2 _ (Al t q-^MC(1 + If)) 2 (^JL.q-^M C (l~|)) 2 , _±_ m 



taL) 2 (vftatY 



i .sys 



where A is respectively either up-going or vertical and B is down-going or horizontal. For the UTM sample, 
the calculated flux is used in place of the number of events. The parameter denoting normalization is a*, 
and for each ratio we introduce a systematic uncertainty parameter e^, weighted by the estimated size of 
the uncertainty, <7i, sys . 

The three xf are summed to form a total Xtot w i tn three degrees of freedom. The value of \\ot * s 
used as a hypothesis test for the cases of v^ «-> v T and v^ <-> v s oscillation. We exclude regions in the 
sin 2 29 — Am 2 plane at the 90(99)% C.L. if the value of % 2 is greater than 6.3(11.3) for % 2 of three degrees 
of freedom. Figure g shows separately the excluded regions for these three alternative oscillation modes, 
along with the allowed region from FC single ring event analysis [J One sees that the parameters allowed 
by the FC data in the sin 2 29 — Am 2 plane are excluded at the 99% confidence level by the independent 
tests for both positive and negative Am 2 sterile neutrino oscillations. 

IV. SUMMARY AND CONCLUSION 

In summary, we have presented three independent data samples that discriminate between the oscil- 
lations to cither tau neutrinos or sterile neutrinos in the region of mixing angle and Am 2 preferred by 



*Although stopping upward-going muons and PC events have similar parent neutrino energies, the limited zenith 
range of — 1 < cos S < with a relatively large uncertainty at the horizon B limits their usefulness and they are 
not considered. 

^The two cases of Am 2 > and Am 2 < are treated continuously Jl2|, using the minimum \ 2 f° r the v s fit 
which was slightly lower for the case of Am 2 > 0, and this was used to draw the contours for both cases. 



the majority of the Super-Kamiokande data. Two-flavor oscillation between muon neutrinos and sterile 
neutrinos fit the low energy charged current data, but do not fit the neutral current or high energy data. 
We cannot exclude more complicated scenarios in which both v^ L <-> v T and v^ «-> v s oscillations co-exist 
with small mixing to sterile neutrinos, or with much smaller mass difference for sterile neutrinos; yet 
there is nothing in this data to encourage one about the existence of sterile neutrinos. Pure v^ «-> v r 
neutrino oscillations fit all of the data samples presented, without any inconsistency. 
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FIG. 1. (a,c,e) Zenith angle distributions of atmospheric neutrino events satisfying cuts described in the text: 
(a) multi-ring sample, (c) partially contained sample, and (e) upward through-going muon sample. The black 
dots indicate the data and statistical errors. The solid line indicates the prediction for v^ <-> v T , and the dashed 
for Vy, <-> u s , with (Am 2 , sin 2 2#) = (3.2 x 10 -3 eV 2 ,l). The two predictions are normalized by a common factor 
so that the number of the observed events and the predicted number of events for Vy, <-► v T are identical. (b,d,f) 
Expected value of the corresponding test ratio as a function of Am 2 . The solid horizontal liness indicates the 
measured value from the Super-Kamiokande data with statistical uncertainty indicated by dashed lines. Black 
dots indicate the prediction for v^ <-+ v T , and empty squares for z/ M <-» i/ s , in both cases for maximal mixing. 
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FIG. 2. Excluded regions for three oscillation modes, (a) v^ <-► v T , v^ <-► v s with (b) Am 2 > and (c) Am 2 < 
0, the light(dark) gray region is excluded at 90(99)% C.L., Thin dotted(solid) line indicates the 90(99)% C.L. 
allowed regions from the analysis of FC single-ring events. 



